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Abstract. In this paper, we use the methods found in 1211 to create a con- 
tinuous variable analogue of Shor's quantum factoring algorithm. By this 
we mean a quantum hidden subgroup algorithm that finds the period P of a 
function 

$ : K > R 

from the reals R to the reals R, where <I> belongs to a very general class of 
functions, called the class of admissible functions. One objective in creating 
this continuous variable quantum algorithm was to make the structure of Shor's 
factoring algorithm more mathematically transparent, and thereby give some 
insight into the inner workings of Shor's original algorithm. This continuous 
quantum algorithm also gives some insight into the inner workings of Hallgren's 
Pell's equation algorithm. 

Two key questions remain unanswered. Is this quantum algorithm more 
efficient than its classical continuous variable counterpart? Is this quantum 
algorithm or some approximation of it implementable? 



Contents 

1. IntroductiorJ 2 

2. Mathematical machinervl 3 
2jl^_^en£rali^e^jun£^i£M 3 

2.2. Rigeod Hilhort sr^a.cej 3 

3. Fourier analysis on the real line M 4 

4. The algorithm for finding integer periodsl 5 

5. The observable d 7 

6. The algorithm for finding rational period j 8 
^^^^jd^ ^^^^^g^p^oda ^ 9 
8. Conclusion! 10 



Date: March 15, 2004. 



2000 Mathematics Subject Classification. Primary 81P68; Secondary 68Q05, 81Q99. 

Key words and phrases. Quantum algorithms, continuous variable, Shor's algorithm, general- 
ized functions, distributions, rigged Hilbert spaces, Gel'fand triplets. 

This effort is partially supported by the Defense Advanced Research Projects Agency (DARPA) 
and Air Force Research Laboratory, Air Force Materiel Command, USAF, under agreement num- 
ber F30602-01-2-0522, the National Institute for Standards and Technology (NIST), and by L-O- 
O-P Fund Grant BECA2002. The U.S. Government is authorized to reproduce and distribute 
reprints for Government purposes notwithstanding any copyright annotations thereon. The views 
and conclusions contained herein are those of the authors and should not be interpreted as neces- 
sarily representing the official policies or endorsements, either expressed or implied, of the Defense 
Advanced Research Projects Agency, the Air Force Research Laboratory, or the U.S. Government. 
The first author gratefully acknowledges the hospitality of the Mathematical Science Research In- 
stitute , Berkeley, California, where some of this work was completed. Both authors would also 
like to thank Howard Brandt and Dan Gottesman, for some helpful conversations, and the referee 
for some helpful suggestions. 

1 



2 



SAMUEL J. LOMONACO, JR. AND LOUIS H. KAUFFMAN 



^^^Aj3 2cndix^_Nunibcr theoretic probabilities. I 10 
10. Rcfcrci'^ 12 



Reference! 12 



1. Introduction 



In this paper, we create a continuous variable analogue of Shor's quantum fac- 
toring algorithm. This algorithm is called a continuous variable Shor algorithm 
for the following reason. Recall that Shor's quantum factoring algorithm j2H|, |27j . 
fF3 reduces the task of factoring an integer N to that of finding the period P of a 
function 

$ : Z — > ZmodiV 

from the integers Z to the integers Z modulo TV. So by a continuous variable 
analogue to Shor's factoring algorithm, we mean a quantum algorithm that finds 
the period P of a function 

$ : K — >R 

from the reals K. to the reals M. 

One of the objectives in creating this continuous variable quantum algorithm was 
to make the structure of Shor's factoring algorithm more mathematically transpar- 
ent, and thereby give some insight into the inner workings of his original quantum 
factoring algorithm |27| , |28| . This continuous quantum algorithm also gives some 
insight into the inner workings of Hallgren's Pell's equation algorithm |17j. 

Whether or not this quantum algorithm is more efficient than its classical con- 
tinuous variable counterpart remains to be determined. By allowing continuous 
variables, the complexity class of problems can easily change. For more insight into 
this issue, we refer the reader to Bartlett et al Moreover, the implcmcntabil- 

ity of this continuous variable quantum algorithm, or an approximation there of, 
also remains to be determined. 

Continuous variable algorithms for two other quantum algorithms are to be 
found in the open literature. A continuous variable analogue of Grover's algo- 
rithm was constructed by Pati, Braunstein, and Lloyd in 0; and a continuous 
variable Deutsch-Jozsa algorithm was recently created by Pati and Braunstein in 

M 

There is also a great deal of literature written on many other areas of continuous 
variable quantum information science. For example, work on continuous variable 
teleportation can be found in [7], IHIj and ^Jj, on continuous variable quantum 
secrecy sharing in [50] and on continuous variable entanglement in j^, and 
on continuous variable quantum error correction in ^ and |15j . 
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To create a continuous variable analogue of Shor's algorithm, we will need to 
make use of the mathematical machinery of generalized functions (also known as 
distributions) and of rigged Hilbert spaces (also known as Gel'fand triplets.) 
For more in depth discussions of this mathematical machinery, we refer the reader 

to m, [El, [12], [H], EHl, and 



2.1. Generalized functions. In regard to generalized functions, the reader is no 
doubt familiar with one generalized function, namely, the Dirac delta function 

S{x) 

on the reals M. We will also make use of the following generalized function 

which is an infinite sum of Dirac delta functions over the lattice {p : n e Z}, where 
P is a nonzero real number. 



2.2. Rigged Hilbert spaces. Wc will make use of the rigged Hilbert space 
with orthonormal basis 

{\x) :xeM}, 

where by orthonormal we mean there is a bracket product on Hk defined by 

{x\y) = S{x-y) . 
The elements of T^r are formal integrals of the form 



dx f (x) \x) , 

where / : R — > C is a function or a generalized function. 
For xq a constant, we define 

oo 

\xo) = I dx 5 [x — xo) \x) 



Since the Dirac delta function is a tempered distribution j25|. it follows that 

1 if xo = yo 



(2/0 1 a;o) 



otherwise 
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3. Fourier analysis on the real line M 



Let $ : K — > M be a periodic admissible function of minimum period P from 
the reals R to the reals M. 



Remark 1. We have intentionally not defined the term, 'admissible,' since there 
are many possible definitions of this term. For example, one workable definition 
of an admissible function is a function that is Lebesgue integrable on every closed 
subinterval of the reals M. 



We seek to define the Fourier transform of <I>. Since $ is in general neither 
nor nor of compact support, the usual definitions of the Fourier transform will 
not apply. So we need to be a bit creative. 



We proceed to define the Fourier transform as follows: 



Definition 1. Let $ : R > R be a periodic admissible function of minimum, 

period P from the reals M to the reals R. We interpret the standard expression 



oo 



dx e 27^*2:!/$ ^x) for the Fourier transfc 



Wm $ : R 



C as the generalized func- 



— oo 



tion 



p 








where 




and where C denotes the complex numbers. 
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Remark 2. The above definition can he motivated as follows: 

oo ^ («+l)P 

J dx e'^''"y^{x) = J2 J dx e-^'''''y<i>{x) 



OO „ 

^ jdx e-2'^^(^+"-P)!/$ (x + nP) 



n — — OOr 



OO „ 





p 



= w]^ (y - 7) e^2."^$ [x] 



= Sp (y) J dx e^^TTixy^ ^^j) 


where, in the context of distributions, we have 

CO 

n— — 00 ' ' 

(^^ee Eg.; 



TO TO + 1 



The reader can easily verify that the inverse Fourier transform behaves as ex- 
pected, i.e., that 

Proposition 1. 

00 

$ (x) = J dy e-^'^"^* (y) 



4. The algorithm for finding integer periods 



Let 

$ : M — >R 

be a periodic admissible function of minimum period P from the reals M. to the 
reals M. We will now create a continuous variable Shor algorithm to find integer 
periods. In later sections, we will extend the algorithm to rational periods, and 
then to irrational periods. 



We construct two quantum registers 
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I Left Register) and | Right Register) 

called left- and right-registers respectively, each 'living' respectively in its own sep- 
arate rigged Hilbcrt space TYr. The left register is constructed to hold arguments 
of the function the right to hold the corresponding function values. 



We assume we are given the unitary transformation 



defined by 



U $ : Tim. <8) — > ® Tii 



U^:\x) \y)^\x) \y + ^{x)) 



Finally, we choose a large positive integer Q, so large that Q > 2P^. 



The quantum part of our algorithm consists of Step through Step 4 
described below: 



as 



Step Initialize 



l^o) = |0) |0) 



Step 1 Apply the inverse Fourier transform to the left register, i.e. apply ^ ® 1 
to obtain 



IVi) = / dx e^'^"-'^ \x) |0) = / dx \x) |0) 



Step 2 Apply : \x) \u) i — > \x)\u + ^{x)) to obtain 



1^-2) 



OO 

J dx \x) 
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Step 3 Apply the Fourier transform to the left register, i.e. apply (g) 1 to obtain 



oo oc 



(n+l)P 



IV'3) = dy dx 6-2--!/ \y) \<p{x)) = dy \y) J2 / <r^^"'y \^ {x)) 



-oo — oo 



nP 



[ dy \y) /t^a; e-2'^^(^+"-P)^ 1$ (x + nP)) 

-oo "=-°° 

Jdy \y) ( ^ e-^""^^"] ( J dx e'^"-^ |$ {x)) 

-oo \n=-oo / \o ^ 

Jdy\y)Sp{y)(jdxe-'^'^y\^{x))\= £ | ^) f ^ Jcix e'^'^^-f |$ (x)> 



where 



n 



Step 4 Measure the left register with respect to the observable 

oo 

0= jdy^\y){y\ 



to produce a random eigenvalue 



m 

Q 



where \_Qy\ denotes the greatest integer < Qy, and then determine whether 
^ can be used to find the period P. 



5. The observable O 



In this section, we now discuss the above Step 4 in greater detail. 
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The spectral decomposition of the observable O is given by 

oo 

IQx 



O 



dx 



m— — oo 



where Pm denotes the projection operator 



Pm = j dx \x) {x\ 



Measurement of the left register of 

oo 

IV'3)= E 



n— — oo 



p 



with respect to O will always produce an eigenvalue ^ for which there exists an 
integer n such that 

m n m+1 

— < — < 

Q - P Q 

We seek to determine the unknown rational from the known rational eigenvalue 
Q- 

If Q > 2P^, then the unknown rational ^ will be a convergent of the contin- 
ued fraction expansion of the known eigenvalue ^. Thus, the continued fraction 
recursion can be used to determine the period P. (See jlBI, Theorem 184, Section 
10.15].) 



6. The algorithm for finding rational periods 



We now extend the above algorithm to one for finding rational periods 



P=-, gcd(a,6) = l 





We choose an integer Q > 2a^ 



Part 1 Execute the above steps Step through Step 4 twice to produce two 
eigenvalues 

mi 7712 



and then goto Part 2 



Since Q > 2a^, the eigenvalues ^ and ^ will have unique convergents respec- 
tively of the form 

nib 7126 
and 



(See [H Theorem 184, Section 10.15].) 
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If the following Condition A is satisfied, then the reciprocal period is simply 
given by 

1 gcd(ni6, n26) 



Condition A. gcd(ni,n2) = 1, gcd(ni,a) = l, gcd(n2,a) = l 



If we assume that Condition A is satisfied, then the above expression for the 



reciprocal period can be computed in Part 2 given below: 



Part 2 



Step 5 



Execute the following: 

Compute all the convergents : = 1, 2, 



Step 6 



and respectively 



, and | 



^-.£=1,2, 



Search for denominators qik and q2e which are equal 
For k = 1,2, . . . ,K do 
For£ = 1,2,...,L do 
If qik = q2i then 

Let q = git. = 02/ and a = — ^r-^ r 

If a is a period of $ Then 



Output a and Stop # Period found 



EndFor 
EndFor 

GOTO 



Part 1 # Period not found 



Part 2 will find and output the period P provided the output of Part 1 



satisfies Condition A. 
will occur after 



Part 1 



From the last corollary of the Appendix, wc know this 
is repeated an average of O ((Iglga)^ ) = O ((IglgQ)^ 



times. However, since we do not know until the completion of Part 2 whether 



or not the output of Part 1 satisfies Condition A , both Part 1 and 



Part 2 



need to be repeated on average at most O ^(IglgQ)^^ to finally find the output P. 



Remark 3. One can quadratically speedup Step 6 by taking advantage of the fact 



that the convergent denominators are linearly ordered. 



7. Finding irrational periods 



The above algorithm can be extended to finding, to any degree of desired pre- 
cision, the period P of a periodic admissible function when the period P is 
irrational. But in this case, there is a severe restrictive condition that must be 
imposed on the function Namely, we need to assume that the function $ is 
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continuous. This continuity condition is needed for determining whether or not a 
rational is sufficiently close to the unknown irrational period. 



8. Conclusion 



The continuous variable quantum algorithm constructed in this paper does give 
some insight into the inner workings of Shor's original quantum factoring algorithm. 
Moreover, it also gives some insight into the inner workings of Hallgren's Pell's 
equation algorithm 17). 

On the other hand, the quantum algorithm constructed in this paper raises many 
more questions than it answers. Is this quantum algorithm more efficient than its 
classical continuous variable counterpart? Can this algorithm be implemented? 
Can an approximation of this algorithm be implemented? 



9. Appendix. Number theoretic probabilities. 



In this Appendix, we derive an asymptotic lower bound 



Is ig 



the 



probability that the output of 



Part 1 



of the algorithm found in Sectional of this 



paper will satisfy the Condition A defined within that Section. 



Notation Convention. Throughout this section, the symbol 'p ' will always be 
used to denote a prime integer. 



Proposition 2. Let a be a fixed positive integer. Then for every positive integer 
N > a, if an integer n is randomly chosen from the set integers 

{keZ:0<k<N} 

according to the uniform probability distribution, then the probability 



ProbN l^gcd (a, n) = 1 
that n is relatively prime to a is bounded below by 
ProbN (gcd (a, n) = 1 ) > 



a 

where Z denotes the set of integers, and where ip denotes the Euler phi function. 
Proof. 

p™v u („, „) = i) . n (i - > n (i - - " 

p\a p\a 



□ 



As a corollary, we have: 
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Corollary 1. Let a be a fixed positive integer. Then for every positive integer 
N > a, if ni andn2 are two random integers chosen independently with replacement 
from the set integers 

{keZ:0<k<N} 
according to the uniform probability distribution, then the probability 



ProbN ygcd (a, ni) = 1 = gcd (a, 
that both Hi and n2 are relatively prime to a is bounded below by 

ProbN I gcd (a, tii) = 1 = gcd (a, 712) ) > 



a 

where TL denotes the set of integers, and where ip denotes the Euler phi function. 

Proposition 3. Let a be a fixed positive integer. Then for every positive integer 
N > a, if ni andn2 are two random integers chosen independently with replacement 
from the set of integers 

{keZ:0<k<N} 
according to the uniform probability distribution, then the conditional probability 



ProbN |^gcd(ni,n2) = 1 | gcd (a, ni) = 1 = gcd(a, ^2) 

that ui and 712 are relatively prime given that ni and n2 are both relatively prime 
to a is bounded below by 

ProbN ^gcd (ni, 712) = 1 | gcd (a, m) = 1 = gcd (a, 712)^ > 

where Z denotes the set of integers, and where ip denotes the Euler phi function. 
Proof. 



ProbN gcd (711,712) = 1 



ycd (a, 712) = 1 / p 

p\a and p<N 



V N 



> n i'-p-') 

p 

p'fa and p<N 

where C denotes the Riemann zeta function. (See ^Hl-) ^ 



Corollary 2. Let a be a fixed positive integer. Then for every positive integer 
N > a, if ni and 712 are two random integers chosen independently with replacement 
from the set of integers 

{keZ:0<k<N} 
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according to the uniform probability distribution, then the probability 

ProbN ^gcd(ni,n2) = gcd(a,ni) = gcd(a,n2) = 1 
that the integers a, ni, n2 are all relatively prime to each other is bounded below by 

ProbN ^gcd (rti, 7^2) = gcd (a, ni) = gcd (a, 712) = > 

where Z denotes the set of integers, and where tp denotes the Euler phi function. 
Moreover, we have the asymptotic bound 

ProbN [ gcd (711,712) = gcd (a, ?ii) = gcd (a, 712) = 1 ) = 



6 f Lp(a)^ 



Iglga 

Proof. The first part of this corollary follows immediately from the above corollary 
and proposition. The second part follows immediately from a number theoretic 
theorem found in [181 Theorem 328, Section 18.4] which states that 

liniinf^i^ = e-^ 
a/ In ma 

where 7 denotes Euler's constant. □ 
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